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Abstract. In an attempt to characterize the distribution of forms and shapes of nodal 
domains in wave functions, we define a geometric parameter - the ratio p between the 
area of a domain and its perimeter, measured in units of the wavelength \j\f~E. We 
show that the distribution function P{p) can distinguish between domains in which 
the classical dynamics is regular or chaotic. For separable surfaces, we compute the 
limiting distribution, and show that it is supported by an interval, which is independent 
of the properties of the surface. In systems which are chaotic, or in random-waves, 
the area-to-perimeter distribution has substantially different features which we study 
numerically. We compare the features of the distribution for chaotic wave functions 
with the predictions of the percolation model to find agreement, but only for nodal 
domains which are big with respect to the wavelength scale. This work is also closely 
related to, and provides a new point of view on isoperimetric inequalities. 



1. Introduction 

In this work we study the (real) eigenf unctions of the Laplace-Beltrami operator — Am 
on a Riemannian surface M. with Dirichlet boundary conditions (if A4 has boundaries). 
Consider a real eigenfunction which satisfies 

(A M + E J )4j J (r) = , i/> j (r)\ r edM = Q - (1) 
The nodal domains are the maximally connected domains in M. where ipj has a constant 
sign. The nodal set (or the set of nodal lines) is the zero set: Uj — {r G M. : ipj(r) = 0}, 
which also forms the boundaries of the nodal domains. We shall denote the nodal count 
(i.e. the number of nodal domains) of if>j(r) by Uj. 

The investigation of quantum signatures of classical chaos and integrability has been a 
hot topic in quantum chaos for a long time [TJ [2]. In the past few years, the interest 
in nodal domains, their counting and their morphology increased after Blum et al [3] 
proposed a quantitative method which distinguishes between the distributions of nodal 
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counts in domains where the underlying classical dynamics is integrable (separable) or 
chaotic. This added a new approach, the statistical investigation of nodal patterns, to 
the more common investigation methods of spectral or wavefunction statistics, which are 
often connected to random-matrix theory [I]. Blum et al showed that if Uj is the nodal 
count of the jth energy eigenstate of a domain , then £j- = Vj/j has a limiting distribution 
-P(£) = lim n ^oo - Y^j=i 5 where the characteristics of the distribution depend on 

the classical properties of the domain. For separable domains, P(£) has a square root 
singularity at a (system-dependent) maximum value, while for chaotic systems P(£) 
is (approximately) normally distributed. Comparison between the numerical results 
for chaotic billiards and the random-wave ensemble supports Berry's conjecture [5] - 
wave functions in a chaotic system behave in the limit of high energy like a random 
superposition of plane waves. By that, the qualitative observation of Miller et al (6J, that 
nodal sets can be used to distinguish between wave functions in chaotic and integrable 
domains, could be tested in a quantitative way. Other studies of various quantities - 
which pertain to the morphology and complexity of the nodal network - were published 
in the mathematical and physical literature, building upon the older results regarding 
the bounds on the total lengths of the nodal lines and their curvature E]. E.g., in [9] 
the distribution of the curvature is calculated, in addition to the mean and the variance 
of the total length of the nodal set. The distribution of the avoidance distances between 
nodal lines was also computed [10J, to mention few examples. 

An important breakthrough has been achieved by Bogomolny and Schmit [11] who 
implemented a critical percolation model that explains the large scale structure of 
nodal domains in chaotic wave functions. This model is supported by a variety of 
numerical calculations. For example: the expectation value and variance of the nodal 
count for chaotic billiards, as well as the distribution of areas of nodal domains, follow 
the predictions of the model [31 [11]; The nodal lines in the high energy limit seem 
(on large scales) to be SLE 6 curves p21 021 EE] as it is proved for the boundaries of 
percolation clusters [15]. Despite the good agreement, the percolation description is a 
priori insensitive to the structure of the nodal set in scales of the order of a wavelength. 
In addition, it was demonstrated by Foltin et al [16] that there are some special measures 
with a scaling behavior which is different for percolation and the nodal set of the random- 
wave ensemble. The latter special measures, in general, probe subwavelength scales at 
two points at a (large) distance. 

In this work we suggest a new (quantum mechanical) method for the classification 
of billiards according to their classical properties. We will discuss below in what sense 
the signatures in nodal patterns differ from the scenario known for the more common 
spectral and wavefunction analysis. 

Our method provides yet another test to the conjectures by Berry and Bogomolny. The 
parameter which we use in order to interrogate the morphology of nodal lines is defined 
as follows - We consider the jth eigenfunction of and its nodal domains sequence 
{tjJj} m =\... Vj - The indices j, m specify a nodal domain; For this domain we define the 
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r»' ~~ r (m) • y A ) 

where A^ and Lj are the area and perimeter of and the ratio is measured in 
units of the wavelength 1 / y/E. We shall define for different ensembles two different 
probability measures on the parameter p. 

For wave functions which satisfy JT| on a compact domain A4, we consider a spectral 
interval I = [E,E + gE], g>0 with JV> = (J{j : £y G 1} , and define: 

PM( P ,E,g) = ±J2^Y, 5 (p-Pj n) )- ( 3 ) 

J J m=l 

Note that in the above, the weights of nodal domains which belong to the same 
eigenfunction are equal, but not necessarily the same as the weight of domains which 
belong to another eigenfunction. 

The second probability measure pertains to an ensemble of wave functions on unbounded 
domains, in our case - the gaussian random-wave ensemble (which will be described in 
section [3]). Since the wave functions do not satisfy any boundary condition, we consider 
them over an arbitrarily large and fixed domain Q C K 2 , and include only the nodal 
domains which are strictly inside Q. We denote their number for a given member of the 
ensemble by uq, and define: 



P rw (p,E,n) = ( — Y6(p-p 3 )) (4) 



The reason for using two different measures is the different nature of the problems at 
hand. However in the limit the two measures coalesce. We shall investigate the existence 
and the features of a high energy limiting distribution 

P(p)= lim P(p,E). (5) 

E— >oo 

The choice of the area-to-perimeter ratio p as a parameter to characterize the geometry 
of nodal domains is is inspired by the following considerations: The nodal pattern 
for separable surfaces is a checker-board, where a nodal domain uj is asymptotically a 
rectangle with sides of the order of a wavelength. Therefore ~ E~ l ,L U ~ E~^ 
and Pu will be of the order of one. Similarly, according to the percolation model, a 
nodal domain of a chaotic surface is asymptotically shaped as a chain [11] with n cells, 
where for each cell A c ~ -E -1 , L c ~ E~% where L c is the cell's contribution to the nodal 
domain's perimeter (see fTTI) . We get that in both cases the parameter p for a typical 
nodal domain will be of the order of unity, yielding localized distributions for the two 
types of surfaces. However, as will be shown below, the distributions differ substantially 
for systems with different classical properties. 

In addition, the area-to-perimeter ratio p is relevant not only to the study of the high 
energy limit, but arises as a natural parameter in the study of isoperimetric inequalities 
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(see e.g [17J). The restriction of a wave function ipj to one of its nodal domains uo. 



m ■ 



is an eigenfunction of the Laplace-Beltrami operator on the domain u m , with Dirichlet 
boundary conditions. Since it consists of a single nodal domain, Courant theorem [18] 
implies that it is the ground-state of u m . Therefore, knowing pj we can express the 
ground-state energy in terms of the area and perimeter of u m . In the mathematical 
literature there are known bounds for such expressions - A relevant example is the 
bound for convex domains, derived by Makai [19] (lower bound) and Polya [20J (upper 
bound, which was generalised to all simply or doubly connected domains by Osserman 



In order to derive a distribution function P(p) between the extreme values, some measure 
on domains should be defined. Here we confine ourselves to well defined families 
of domains - those obtained as nodal domains of a given ensemble, and due to this 
restriction we are able to define the measures ([3]) , (j4j) and study the limiting distribution 
for different classes of systems. 

In this paper we shall examine the distribution function P(p) for separable and chaotic 
domains and for the gaussian random-wave ensemble. We will show that 

• The limit distributions we obtain have strong "universal" features. That is, they 
depend crucially on the type of classical dynamics the manifold supports, and only 
to a lesser extent on the idiosyncratic details of the actual system. 

• The limiting distribution for the random-wave ensemble is similar to the one for 
chaotic domains, as predicted by Berry's conjecture. 

• The limiting distribution for random-waves (chaotic billiards) is consistent with the 
percolation model but contains (universal) information beyond percolation as short 
length scales on the order of a wavelength are probed for small nodal domains. 

The limitation to (quantum mechanically) separable systems is due to the checkerboard 
structure in the nodal patterns of their wavef unctions. A generalization to all integrable 
or pseudo-integrable domains, where, in general, the checkerboard structure is lost, 
would be desirable. This highlights a general difference in the scenario known e.g. from 
spectral statistics where all integrable systems (separable or not) share the same 
(Poissonian) statistics. Quite contrary statistical properties of nodal domains in 
integrable systems are very different for seperable systems with a checkerboard structure 
and non-separable systems with generically no nodal crossings pTJ|, [22] . Note, that 
when a separable system is slightly perturbed, all nodal crossings will open (in an often 
highly correlated way which makes the introduction of a percolation model at this point 
quite difficult) and statistical properties will change singularly (again in contrast to 
what is known from spectral statistics where a small perturbation smoothly changes the 
statistics). 



EH): 



7T 7T 



(6) 
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2. The limiting distribution of the area-to-perimeter ratio for separable 
domains 

As was mentioned above, the nodal network of eigenf unctions of separable surfaces has a 
checkerboard structure. This follows from the fact that one can always choose a basis in 
which all the eigenfunctions can be brought into a product form. Therefore, one might 
expect that the main features of the limiting distributions of different surfaces will be 
similar. We will show that this is the case, and therefore we begin this section by explicit 
calculation of P re c{p) for a rectangular billiard. The discussion of this simple example 
will pave the way to computing P(p) for other systems (the disc billiard and a family 
of surfaces of revolution) and to the identification of some common features which we 
assume to be universal for all the separable systems. The detailed computations are 



presented in Appendix A 



The Dirichlet eigenfunctions for a rectangular billiard with side lengthes a, b are: 

4>mn{x, y) = sin sin — — = ipm{x)ipn{y) (7) 

a b 

The corresponding eigenvalues are: 



E = TT 



fm\ 2 (n^ 2 



E m + E n (8) 



where E m , E n can be interpreted classically as the energy stored in each degree of 



freedom (A formal definition can be found in Appendix A ). The nodal domains are 



rectangles of size: -7= x -7==, therefore: 

V J^m V 



7T 2 _ 1 I 

Ann — z=r- i L mn 2,71 ( __ — -\- 



and 



-1 



Pmn = l /^L + Jl-*=-) 3 ,(*S.) (9) 
2 \ V E mn V E mn I y E mn J 

where: p fjr^H = P (ir 22- )- ^ ne mere form of ([HD implies that p is bounded by: 

-j= ^ p ( Z c) ^ \ (10) 

Thus, the support of the distribution function P(p) is an interval which is narrower than 
©■ 

In the limit E — > 00, the distribution ([3]) can be approximated (neglecting corrections 
of order E~*) by an integral. Performing the integration over the variables: 

En 



k = \/E mn , 9 = arctan(— ^ 

E, 



in 



the distribution function is: 



2 r* ( 7T \ f 7 4 for -is ^ p ^ f 

p rec ( P ) = - Hp- * = { pV ^ n . ( n ) 

V 2(sm6 + cos9)J Q otherwise 



7T 
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Figure 1. The limiting distribution for a rectangular billiard (fTT|) compared to the 
calculated distribution for eigenfunction with E ■ A < 10 5 



The explicit form of P rec (p) suggests the following qualitative and quantitative 
conclusions: 

(i) The existence of a limiting distribution function (which is independent of the aspect 
ratio a/6 of the billiard) is demonstrated. 

(ii) It is supported by the compact interval [7r/\/8, 7r/2] . 

(iii) P re c{p) is a n analytic and monotonic decreasing function in the interval where it is 
supported. 

(iv) Pr ec (7r/x/8 + 5)^ o+ ~ 1/VS. 

(v) P rec (n/2) = 8/n 2 , Hence, P rec (p) is discontinuous at both boundaries of the support. 

The fact that p mn depends solely on the partition of the energy between the modes 
was a key element in the construction above, and plays a similar role in computing 
P(p) for the other separable systems. When p = we get from equation ([9]) that 

E m = E n , while p = it/2 means that the energy is concentrated completely in one degree 
of freedom. The concentration of probability near p = n/\/8 shows that equal partition 
of energy is prevalent among the nodal domains. 

An explicit derivation of the limiting distribution P(p) for the family of simple 



surfaces of revolution (following Bleher [23]) can be found in Appendix A in addition 
to a separate derivation for the disc billiard. In both cases it is proven that: 

Psep(p) = Prec(p) ' T( P ) (12) 

Where P rec is given by (fTTj) . and T(p) is a finite, positive and smooth function of p. 
Therefore the features which characterize P rec (p) dominate P sep (p) for all the systems 
considered; Thus P sep (p) is supported on the same interval and demonstrates the same 



Geometric characterization of nodal domains: the area-to-perimeter ratio 



7 



type of discontinuities at its boundaries. 

Following the striking similarity of the distributions for all of the investigated manifolds, 
we suggest that properties (i)-(v) of P re c(p) which were derived for the rectangular 
billiard, are universal features of P sep (p) for all two-dimensional separable surfaces. We 



support this assumption by a heuristic model which is presented in Appendix A 




Figure 2. Semi-Classical approximation of the limiting distribution for the disc 
billiard compared with numerical results for 3 different energy intervals: 200 < kR < 
210, 400 < kR < 410 and 600 < kR < 610 where k = \[E 



Numerical simulations for the rectangle and the disc billiards for several energy intervals 
show good agreement with the analytic derivation (see figures fTH2l) . For numerically 
obtained results at finite energies, two kinds of deviations from the limiting distributions 
can be observed: 

(i) Fluctuations along the entire range of p (for the disc) or discrete jumps (for the 
rectangle) in the value of P(p), which vanish in the limiting distribution due to the 
convergence of the corrections to the semi-classical approximations (i.e. turning 
sums over quantum numbers into integrals and the neglect of terms of order E~^). 

(ii) Cusps near p = | and p = | for the disc: the origin for the appearance of these 
features is due to nodal domains with exceptional geometry - the inner domains of 
all wavefunctions (which are asymptotically triangles) for the former, the domains 
of \I/ n .o (which are ring shaped) for the latter. 

It is verified (analytically and numerically) that these differences converge to zero as 
E~z, or faster. 
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3. The limiting distribution of the area-to-perimeter ratio for the 
random-wave ensemble and chaotic domains 

While for chaotic wavefunctions there is no known analytic expression for the nodal 
lines, we will use known results about the morphology of the nodal set in order to 
propose some physical arguments for the expected distribution. The explanations we 
propose are all in agreement with numerical simulations - a detailed information about 
the numerical techniques and the reliability of the results can be found in appendix B. 
A frequently used model for eigenfunctions in a chaotic billiard is that of the Gaussian 
random-wave ensemble. This is based on a conjecture by Berry [5] that eigenfunctions 
of a chaotic billiard in the limit of high energies have the same statistical properties as 
the Gaussian random-wave ensemble. 

A solution ip for the Helmholtz equation (CQ) with a given energy E = k 2 on a given 
domain, can be written as a superposition of functions {^;(r)}^_ 00 which span a 
complete basis, for example: 



Since the solutions of (CD) are real, we are restricted (for this choice of basis) by 
a_i = (— 1) a*. According to Berry's conjecture, expanding the eigenfunctions of chaotic 
billiards (in the high energy limit) in terms of ({TBI , the coefficients a/ distribute for / > 
as independent gaussian random variables with 



and therefore can be modeled statistically by this ensemble of independently distributed 
Gaussian random-waves. 

As suggested by Bogomolny and Schmit [UJ, the nodal domains of a random wave 
are shaped as critical percolation clusters (see Figure [3]), where each site is of an average 
area 



Where (for 2d percolation): r = 187/91. For bond-percolation model over a lattice (as 
illustrated in fig. [3]), the area of a cluster which spreads on n sites is (a\ + a2)n — a2 
where a,\ is the area of a single site and a2 is the area of the connection between two 
sites; the average perimeter is (l± + Izjn — I2 where h, I2 are the average contributions 
to the perimeter of a site and a connection (Since a cluster may contain loops which 
affect its perimeter, we must speak about average). The average area-to-perimeter ratio 



oc 




(13) 



(aia* v ) = 5i >v 



(14) 
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Figure 3. Above - a realization of the bond-percolation model by Bogomolny and 
Schmit. The dark sites are positive, and the bright are negative; at each junction there 
is a saddle which can be positive or negative with equal probability, connecting by 
that two of its neighbors to the same cluster. Below - a typical realization of the nodal 
structure of a random-wave on the wavelength scale. 



can be written as 



L~ Cl 



1 



Co 



(17) 



n + C 3/ 

where C%, C 2 , C 3 > 0. This relation can be used as a guideline to the desired distribution 
of p. 

Indeed, as was confirmed numerically, the distribution P(p) for nodal domains follows 
(j!7p in several aspects. We have examined the restricted distribution for nodal domains 
with a given number of sites - we define P^ (p) to be the distribution for nodal domains 
with area 



n 



l)As < A {n) < 



1 



n 



2' 2' 
We found out that P {n \p) is roughly symmetric about a mean value: (p n ). As in (Q3 
(p n ) is increasing with n, and converging to a limiting value: poo. Since the percolation 
model is assumed to provide an exact description of the system in the high-energy limit, 
we expect ffTTj) to serve as a good approximation for large domains (see fig. [5]). 
The value of p^ is a direct result of a theorem by Cauchy for the average chord length 
of a domain: 

71 ■ A 



(19) 
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Where a is the chord length. The original theorem (which was stated for convex 
domains) is extended in [25], to include nonconvex and multiply connected domains. 
For nodal domains of infinite size, the statistics of the average chord length should follow 
that of the entire nodal set. That in turn is known to be [5l [26]: {(Jrw) — V^n/k, 
therefore: 

Poo = = V2 (20) 

7T 

The value of (pi) can also be estimated: as shown in [10] . the single cell nodal domains 
are mild deformations of a circle of radius r = jo/k (where jo ~ 2.405 is the first zero 
of Jo). Therefore: 

(Pi) = — = J ~ (21) 
Kn/ 2nr 2 V ; 




1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 1.55 
P 



Figure 4. The distribution function P^ n \p) is plotted for several values of n. The 
mean value is increasing with n, while the variance decreases (except for n = 1) as 
suggested by the heuristic model. 

In order to study the impact of the deformations on the value of (pi), we have calculated 

p for a variety of domains, like ellipses, rounded shapes with corners (a quarter of a circle 

or a stadium etc.) and others. The results show that stretching of the nodal domain 

(e.g. increasing the eccentricity of an ellipse) increases p, while turning it "polygonal" 

(i.e. having points on the nodal line with very high curvature) reduces p. 

Derivation of (p n ) for other values of n seems to be more complicated. However, fitting 

between the numerical results for high n values and ( TT7l) equips us with the empirical 

result (which is valid for n ^> 1): 

. . r n + 0.805 , . 

(p n ) w V2 22 

n + 0.936 V ; 
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n* ** * 



^ numerical results 
best fit for high n's 



15 20 25 30 

n 

Figure 5. The measured value of (p n ) compared to the best numerical fitting (for 
high values of n) to (JT7J): (p n ) = y/2(n + 0.805)/(n + 0.936). 

Another interesting feature is the width of the distribution around (p n ). Equation <HM 
implies that the variance is proportional to the variance in the average chord length 
between different nodal domains of the (approximately) same area. Therefore, the 
variance is expected to be smaller for larger domains, which follows the statistics of the 
entire nodal network to a larger extent. The only exception is the variance for single 
site domains, which as was mentioned [10], have strong limitations on their shape, and 
therefore a relatively small variation in the average chord length. 

The bounds ([6]) on p should not hold in general for the nodal domains of (fT3l) ; however, 
the numerical bounds seem to agree with ([6]) for all of the measured nodal domains, 
including multiply connected domains, suggesting that ([6]) is valid for the nodal domains 
of the ensemble with probability 1. 

The distribution of p for all of the nodal domains is given by: 

oo 

P(p) = 5>HP^(p) (23) 

n=l 

where p(n) is given asymptotically by (TT6j) . 

Fig. [6] shows the calculated distribution for 3 different systems - A random-wave 
ensemble, the inner domains of a Sinai billiard and those of a stadium billiard. 
Comparing the functions, we find additional strengthening to Berry's conjecture. There 
are boundary effects of chaotic billiards - e.g. a peak in the distribution P(p) near tt/ v^8 
and \/2, and a lower probability for p ~ jo/2, however they vanish in the semi-classical 
limit (see fig. [7]). In our study it was easy to put those effects aside - if we consider only 
inner nodal domains for the chaotic billiards (as in fig. [6]), we observe no prominent 
differences between the distributions. 
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Figure 6. A comparison between the distribution function P(p) calculated for the 
random-wave ensemble and for inner domains of a Sinai and stadium billiards. 




Figure 7. The distribution function P(p) for a Sinai and stadium billiards. We 
compare the distribution for relatively low energy intervals (asterisks) , higher intervals 
(crosses) and inner domains only (solid line). The distribution function indeed (slowly) 
converges to the distribution of inner domains (the peaks at n/y/8 and y/2 are lowered, 
the one at jo/2 is elevated), which is similar to the distribution for the random- wave 
ensemble (see fig. 11) 



4. Conclusions 

The main results of this work can be summarized as follows: 

(i) The distribution function P(p) of the area-to-perimeter ratio p, distinguishes 
between billiards with separable or chaotic classical limit (see fig. [8]). 

(ii) The distribution (fT2"l) for the examined separable billiards has some universal 
features, such as a common support and a square-root divergence at the lower 
support. In all studied cases the distribution is a mild deformation of the 
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Figure 8. A comparison between the limiting distribution (fTT|) derived for the 
rectangle billiard and the random-wave ensemble. 

distribution that we found for a rectangle. 

(iii) In accordance with the random-wave conjecture, we find numerically that chaotic 
billiards (stadium and Sinai) have a universal limiting distribution P(p), and it 
converges to the distribution found for the random-wave ensemble. By considering 
only the inner nodal domains for billiards, the agreement can be shown also for 
finite energies. 

(iv) The numerical results suggest that for nodal domains of a random wave or of 
eigenfunctions of chaotic billiards, the area-to-perimeter ratio is bounded by (jSJ), 
i.e. 7r/4 < p rw < 7r/2, including nonconvex and multiply connected nodal domains 
(for which these bounds have not been proven), with probability one. 

(v) We examined the percolation model for the nodal set of random waves from the 
perspective of the area-to-perimeter ratio. It is shown that on the wavelength scale 
the geometry of the nodal domains can only be poorly characterized by percolation 
arguments. However, for large domains the geometry can be described by heuristic 
expressions like ( JT71) . which are consistent with percolation theory. 
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Appendix A. Derivation of the area-to-perimeter distribution for some 
separable surfaces 

In this appendix we suggest a heuristic model for the universal features of the limiting 
distributions of the area-to-perimeter ratio for two-dimensional separable domains. The 
model is supported by an explicit derivation of the limiting distribution for the disc 
billiard and for simple surfaces of revolution. 



Appendix A. 1. Universal features of the distribution 

We begin by considering the classical geodesic flow in a two-dimensional compact domain 
(e.g. a billiard). For a separable domain, a trajectory can be specified by its action- 
variables: 

m — Jp\dq\ , n = j)p 2 dq 2 (A.l) 

where q\, q 2 are the (separable) coordinates, pi,p 2 are the conjugated momenta and the 
integration is over one period of the specified coordinate. At every point along the 

1 I • 1 2 

trajectory, the energy E mn = | |q| can be expressed as E mn = E m + E n , where: 

E m = ^\q-qi\ 2 , E n = ^\q-q 2 \ 2 (A.2) 

where qi, q 2 are the local unit vectors - if we consider circular domain for example, then 
qi = r = (cos#, sin#), q 2 = 9 = (— sin 9 , cos 9) . In general E m ,E n are not constants of 
motion. The only exceptions are the trajectories in a rectangular billiard. 

From a quantum point of view, the eigenstates of Schrodinger equation (CQ) for a 
separable domain, can be written as ip mn = VVn^OVvfe) , while the Laplace-Beltrami 
operator can be written as: A = A m + A n , where A m i/)„ = A n ip m = 0. This allow us to 
define the quantum analogue to (1A.2I) : 

j? ( \ A m ^(gi,g 2 ) A n ip(q u q 2 ) 

E m (q u q 2 ) = — , E n (qi,q 2 ) = — (A.3) 

In the semi-classical limit, the spectrum of a separable domain is given by 
{E mn \m,n G N}, where E mn is the energy of the classical trajectory specified by the 
action- variables mh~, nh~ (as emerging from Bohr-Sommerfeld quantization) . In addition 
the semi-classical value of (1A.3j) converges to the classical value (1A.2I) for every point in 
the domain. 

In section [2], the value of the area-to-perimeter ratio for a given realization i/j mn for the 
rectangular billiard, was derived to be 

Pmn = i ( njj?L + a (a.4) 

2 \\ E mn V E mn J 

Therefore, for a rectangular billiard, the value of the (quantum) parameter p, has also 
an immediate classical interpretation (see figures |AT|A2I) . 
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Figure Al. The value of p as a function of the partition of energy between the two 
degrees of freedom, for the rectangular billiard. 




0.2 0.4 0.6 0.8 1 

E„, 
£„,„ 



Figure A2. The limiting distribution of the partition of energy between the two 
degrees of freedom for a rectangular billiard. 

In order to generalize the limiting distribution which was derived for a rectangular 
domain to other separable domains, we suggest the following heuristic model: 

• In the high energy limit, (almost all of) the nodal domains of a separable domain 
are converging to rectangles, and the wave function in the close neighborhood of a 
nodal domain is converging to (J7J). Therefore (in the limit) equation (1A.4I) should 
hold. However, since for general separable domains Vw( r ) is n °t an eigenfunction 
of the operators A m , A n , the value of pmn for a given realization ip mn will depend 
on the interrogated nodal domain ujj . 

• In the continuum limit, the limiting distribution is of the form: 

p{p) = jf-Jgfa n )j M Hp- pSL) ( a - 5 ) 

where the first integral is over the energy interval, and the second is over the domain. 
The function g(m,n) is the quotient of the appropriate Jacobian and v mn . Since in 
the vicinity of p = ir/yS two solutions for (J9j) coalesce, we expect the square root 
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singularity at P(p — > n/\^8 + ) to be a universal feature. 

• In addition, since p{E m j E mn ) is convex (see equation IA.4I and figure lAljl . the 
distribution function should be monotonically decreasing. 

This supports the assumption that the properties (i)-(iv) for p rec and P rec (p), which 
were derived in section [2] for the rectangle, are universal features of P(p) for all two 
dimensional separable surfaces. Moreover, this model suggest - at least for separa- 
ble domains - that a geometric feature of the nodal pattern, i.e. the area-to-perimeter 
ratio of a given domain, can be deduced directly from the underlying classical dynamics. 

The suggested model is supported by an explicit derivation of the limiting area- 
to-perimeter distribution for several separable domains. In these calculations we 
approximate eigenfunctions and eigenvalues by the WKB method. In addition, we 
approximate sums over quantum numbers (see eq. [3]) by integrals and neglect terms of 
order E~^. The error resulting from these approximations is of the order of E~2 and 
therefore converge to zero in the limit. 

The theme of the derivations is similar. The Hamiltonian TC for these systems is 
homogeneous i.e.: 

H{\m, \n) = \ 2 H(m, n) (A.6) 

This implies that the energy of the state ip mn can be expressed as: 

H(m,n) =m 2 n(l,-) = m 2 -h(-) (A.7) 
V mJ \mJ 

Therefore, integration over the quantum number m becomes trivial. We will also use 
the first term in the Weyl series: Jj{j : Ej < E} = AttE/ A + 0(VE) in order to estimate 
Nj. 

Appendix A. 2. The disc billiard 

Equation (OQ) can be written in polar coordinates as 
/ d 2 Id Id 2 \ 

For the disc billiard the boundary condition are: ip\ r =i = 0. The eigenfunction and 
eigenvalues of ( 1A.8[) are: 

ip mn (r, 9) = cos(m6> + ip)J m {j mn r) , E mn = j mn 2 (A. 9) 

where (p is an arbitrary phase and j mn is the nth zero of J m (r). The nodal domains of 
ipmn will be 2m replicas (or one for m = 0) of a slice containing n domains ; we will 
enumerate them as {oJmn}i=i...n, where Um n is the most inner domain. The area and 
perimeter of w™ are: 
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An implicit semi-classical expression for j mn can be deduced by applying the WKB 
approximation to (IA.8I) : 

n = j m yJln-^dr^ (A.ll) 



2 



where: 



I / / m 2 
C mn = - arctan J — I (A. 12) 



7T 



A - 1 



m 2 



it 2 (n+ (§- C mn ) m) 
Setting z = — , z' = — , and substituting (IA.11I) into flA. 12[) we get: 



-•^ + C„-i (A. 13) 

7T Z 

which implies that C mn depends on z solely. Since A mn varies with n as E~^, we can 
approximate: A mn « A m>n _i. Expressing pm„ in terms of (1A.11IIA.12I) yields: 

ffyjmi jm,i — l ) 



Ji) — 2m fA 14") 

m \Jmi ~ Jm,i—1) ~ ^ \Jmi Jm,i—lJ 



7T 



2^ l + sin(27rC(zO) 
Bearing in mind that for a point r G wtL 

= -i-S = + O^- 1 )) (A.15) 



r 2 ip d9 2 j 



in i 



It can be shown that equation (1A.14|) is equivalent to (]A.4|) . 



Integrating (jHJ) over the variables m,C = C(z),C = C'(z'), we get: 



i 



P / (p) = _£ [dmdC-^- f 2 S\p--J 1 , „A dC (A.16) 

W eW/ zdCJc \ H 2 V l + 2sin(27rC")y 1 ; 

Performing the integration we get the limiting distribution: 

P(P) = * = x (A.17) 



n ( 4p 2 + 7Tv/8p 2 -7T 2 sin(27rC) 



where: 



2 \4p 2 - 7rv/8p 2 - 7T 2 7o 1 + (C - |)7rtan(7rC) 

4p 2 -vr v /8p 2 -7r 2 r 2(p) sin(27rC) 
+ 4p 2 + vr v ^ 2 ^ 2= 7o l + (C-i)7rtan(vrC) 

71 ^ = 2^ aiCSin ( " 4p2 4P ) ' 72 ^ = 5 ~ 71 ^ ^ A ' 18 ^ 
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Figure A3. A realization of the functions: Ti(p) = Ci{p) f x .r p \ f( x )dx for the disc. 
T = T\ + T2 is monotonic, finite and positive. 



Appendix A. 3. Surfaces of revolution 

The investigated surface of revolution Ai is generated by a rotation of the analytic 
profile curve f(x) (where x G [—1,1]) around the x axis. We restrict f(x) by: 
f{x) x ^±i = a±\/l -F x i which ensures smoothness at the poles of the surface. In 
addition, we request that f"(x) < 0, so f(x) has a single maximum: f max = f(x max ). 
The Lagrangian of the surface is given by: 

2 



c 



ds 



dt 



(l + f'(x) 2 )x 2 + f(x) 



\2n2 



From which the action variables can be deduced: 
1 



m 



n 



2tt 



f( x ye 



(A.19) 



(A.20) 



2tt 



— (h p x dx = - 



71 



■^(Ef(x) 2 - m 2 )(l + f'(x) 2 )dx 



0. 



where the classical turning points which satisfy: Ef(x±) 2 — m 2 

As in Appendix A. 2 the nodal domains will be 2m + 5 m o copies of a slice with n domains 
and will be denoted by {o>mn}iLi- The homogeneity of the Hamiltonian follows directly 
from (lA.20p . We will follow the notations z = — , z' — — to get: 

(A.21) 



2c2/ 



=^> E(m, z) = m £ 

The WKB approximation to the eigenf unctions is 

cos(m6 + if) 



1pmz(x,6) 



X 



(A.22) 



cos 



4 



+ cos 



kdoc 

4. 
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where 



m 

m 

In the limit of large n, the nodal points density on the curve is high. Therefore, applying 
the WKB approximation, successive nodal points , Xi should satisfy: 

vr = m £ VmW^^^- (A.23) 

In addition, due to the homogeneity of H, x± depend on z solely, therefore: 

iti = m / ^/Elf{xf- 1 V T ^ rfa; => (A.24) 



/ 1 f Xi /mn^i/lEESLj 



z'(xi,z) 



Therefore, the location of the i th zero will be an (implicit) function Xi(z, z') and will not 
be depended on m. The area and perimeter of the nodal domains are: 

J$n = ^ £ mVT+fWdx (A.25) 



£&=^(/C*-i) +/(**)) + P y/1 + f'{x)Hx 



Therefore 



Hm ' n n(f(x t ^) + f( Xi )) + 2m y/1 + f'{xfdx 
Substituting flA~23l) in flA~26l) we get: 

p{z, z>) = I (A.27) 

2l + y /(£(z)f( Xi )y-l 

where £(z)f(xi) = 1 at the turning points. Since for r 6 uomn- 

B ".W = ]$)=7€) (1 + 0( "" l)) (A - 28) 

Equation ( 1A.26I) is equivalent to ( 1A.4I) as well. 
Integrating ([3]) over m, z, z' we get: 

47r Z" 00 r hz) m r 
= —uu\ / rf2 / . / &{p-pWz)W (A-29) 



2tt f 00 cfe 



5(p — p(z ' z))dz 



\M\J Q zS(z)J 

Setting a z / z = £(z)f(x' i ), we get that for p E [ir/y/8, 7r/2] and fixed z, there are two 
allowed values of a: 

2p 



«1,2 



7T 2 — 4p 2 



(tt ± v/8p 2 - 7T 2 ) (A.30) 
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Since f~ l (x) is doubly valued, there are four allowed values of z'\ 



2ty 



z' n = min ( / 1 
for j = 1,2. Therefore 
P(p) 



where 



a 



3 



00 dz 



\M\J Q z£{zf 
P rec (p) ■ (Ti(p) + T 2 (p)) 



E 



z j2 = max ( / 

*V - A) 



-1 



CXj 



p',XA 



-dz' 
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(A.31) 



(A.32) 



T, 



2vr 2 p(4p 2 - n^8p 2 - tt 2 ) 2 



\M\(n 2 - A P 2 ) 2 (tt - y/8p 



71 



27T 2 p(V + 7T^8p 2 - 7T 2 ) 2 




4 2i i=l,2 



|M|(7T 2 -4p 2 ) 2 (7T+ v/8p 2 -7T 2 

£ z is a monotonic increasing function, therefore the integrals in ( IA.32I) diverges at z = 0. 
For 2 >> 1 , = 0(^), therefore the integrals converge at infinity. 
For p — > vr/v^B , the coefficient of the integrals in (1A.32I) has a square root singularity, 
while the integral is converging to a finite positive value, therefore P(ir/ \/8 + 5) ~ 1/Vo. 
For p — > ir/2 — 5, the first coefficient in ( 1A.32I) is of the order of 5. The lower limit of 
integration is defined by: 

df(z') _ S 



Consequently, the first term in (1A.32j) will be of order: 

z'F(z')dz 



6 



(A.33) 



(A.34) 



The value of z' depends on the profile curve f(x), however for z' — > 0, / ~ z'log(z'), 
while for z' — > 00, we get that / ~ l/^' 3 , therefore J is bounded for all possible values 
of z'. 

The value of the second term is due to contributions of nodal domains which satisfy 



Pz 



71 



m 



(A.35) 



2 ■ ■ " f z (z')S 

Therefore, only eigenf unctions for which n/m > 0(1/5) contribute to the second term, 
cLS cL result, it is bounded by: 

J j dndmQin — m5~ l ) Sge 



fj dndm 



(A.36) 



therefore P(tt/2 — 5) is finite, and the universal features specified in section [2] are all 
fulfilled by K32| . 
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Appendix B. Numerical methods for evaluation of the perimeter length 

In order to evaluate the area-to-perimeter ratios and their distribution for the random- 
wave ensemble and chaotic billiards, we have simulated the appropriate wavefunctions 
on a grid. 

We have calculated the statistics for 5000 realizations of random waves, where in each 
realization we summed over 70 terms in (TL3]) ; For chaotic billiards we have reproduced 
the first 2430 eigenf unctions of a Sinai Billiard and the first 2725 eigenf unctions of a 
stadium billiard. 

The (seemingly simple) task of perimeter evaluation must be carried out carefully; it 
can be shown that naive methods, like perimeter's pixels counting, produce an error 
which is independent of the sampling resolution. In order to avoid this error, we have 
approximated the nodal line as a polygon, where the vortices are calculated using a 
linear approximation. We have set the sampling resolution to contain 85 pixels along 
the average distance between two nodal lines (^/2n/k). This resolution was proved to 
produce an error which is less then a percent. The measured perimeter is expected to 
be shorter then the real one, as we are approximating a curve by a polygon. 
The accuracy of the method was tested by calculating the ratio between L^-the total 
length of the nodal set, and the area. It is known that for the random- wave ensemble 

m- 

(L T ) _ 1 

A ■ k 2V2 
The numerical values for this ratio were between: 

0£98 Lr<» W03 

2^2 A-k 2^/2 
It should be noted that when we calculate the total nodal length, we have to calculate 
the perimeter of nodal domains at the edge of the grid. In many cases (i.e. where the 
edge domains are very small) the perimeter calculated for them is larger then the real 
value. It seems likely that the error due to this effect is of the order of the error due 
to polygonal approximation, therefore the two compensate each other, to yield a total 
error which is relatively small, of order ■ 10 -3 . 
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